arXiv: 1502.04704v2 [cond-mat.str-el] 30 Apr 2015 


Memory matrix theory of 
magnetotransport in strange metals 


Andrew Lucas“ and Subir Sachdev“’^ 

“Department of Physics, Harvard University, Cambridge, MA 02138, USA 
^Perimeter Institute for Theoretical Physics, Waterloo, Ontario N2L 2Y5, Canada 


We model strange metals as quantum liquids without quasiparticle excitations, but with 
slow momentum relaxation, and with slow diffusive dynamics of a conserved charge and 
energy. General expressions are obtained for electrical, thermal and thermoelectric trans¬ 
port in the presence of an applied magnetic field using the memory matrix formalism. In 
the appropriate limits, our expressions agree with previous hydrodynamic and holographic 
results. We discuss the relationship of such results to thermoelectric and Hall transport 
measurements in the strange metal phase of the hole-doped cuprates. 


lucasSf as.harvard.edu 
sachdevSg.harvard.edu 


May 1, 2015 


Contents 


1 Introduction 2 

1.1 Transport Without Quasiparticles . 3 

1.2 Outline of the Paper. 4 

2 Hydrodynamics 4 

3 Memory Function Formalism 6 

3.1 Hydrodynamic Interpretation. 7 

4 Dynamics of Charge and Momentum 8 

4.1 Diffusive Transport. 9 

4.2 Momentum Relaxation. 10 

4.3 Magnetic Fields. 12 

5 Dynamics of Charge, Heat and Momentum 13 

6 Thermal and Thermoelectric Transport 15 

7 Conclusions 16 

Acknowledgements 17 


1 









17 


A The Memory Matrix Element 
References 18 


1. Introduction 


The strange metal phase of the hole-doped cuprate superconductors is the most important realization of 
quantum matter not amenable to a quasiparticle description [1, 2]. Apart from the well-known linear in 
temperature (T) resistivity, the strange metal has a Hall angle tan(0H) rsj 1/r^ [3, 4], This combination 
cannot be reproduced in a Boltzmann theory of charge-carrying quasiparticles with a long lifetime; such 
a theory yields a tan(0H) inversely proportional to the resistivity, and so these observations rule out the 
transport of charge by any fermionic quasiparticle, and not just those with the same quantum numbers 
as the electron. 

A general hydrodynamic approach to magnetotransport in strange metals was introduced by Hartnoll 
et al. in [5], and compared to Nernst measurements in the strange metal of the cuprates. The main 
assumption of their theory, apart from the absence of any quasiparticle excitations, was that there was 
a slow mode associated with the decay of the total momentum. Such a slow mode is indeed invariably 
present in proposed field-theoretic models of strange metals [6, 7], and the decay arises from perturbations 
which break the continuous translational symmetry of the field theory (umklapp scattering or impuri¬ 
ties). Combining this momentum mode with reasonable assumptions on the diffusion of the conserved 
U(l) charge and energy densities, general results were obtained by Hartnoll et al. [5] for the charge 
and thermoelectric transport coefficients of a two-dimensional strange metal in the presence of a static 
magnetic field, B. Note that we are using the phrase “charge density” here to refer to the conserved 
density of a global U(l) symmetry. In the application to the cuprates, there are also long-range Coulomb 
interactions associated with such a conserved density, and we will briefly discuss this in Section 4.1. 

Blake and Donos [8] have recently argued that the magnetotransport framework of Hartnoll et al. is 
compatible with the measurements of the longitudinal and Hall conductivities in the cuprates, and also 
provided a solvable holographic model for magnetotransport in the presence of momentum relaxation. 

However, an important concern is that the present derivations of the above results for strange metals 
rely on models which are rather far removed from the microscopic situation in the cuprates. One approach 
[5] begins from a quantum critical point with a relativistic structure, and then breaks Lorentz invariance 
weakly by a chemical potential, the temperature, and the applied magnetic field; the equations were 
then derived using a hydrodynamic gradient expansion restricted by the requirements of the positivity 
of entropy production. The other approach [5, 8] uses holographic models of strange metals represented 
by gravitational theories in a spacetime with an extra dimension, whose field theory duals are not well 
understood. 

In the present paper, we will provide a derivation for the equations for magnetotransport in a normal- 
state metal using the memory matrix approach [9, 10, 11]. Although the memory matrix framework has 
typically been applied to systems where the slow relaxation of momentum dominates magnetotransport 
and the diffusion of conserved quantities is negligible, we generalize and study the case where the diffusive 
dynamics is not negligible and must be consistently included in the memory matrix. In particular, we will 
focus on the consequences of the diffusion of exactly conserved heat, as well as an exactly conserved U(l) 
charge. Our main results may be found in (6) and (70), though there is a large amount of notation that 
must be explained, and so we defer their presentation. Our assumptions on the momentum mode, and 
the charge and energy diffusion will be the same as those in [5], but we will not assume any relativistic 
structure on the underlying theory. Using certain assumptions, our results agree with those of [5], and 
with the recent holographic papers [12, 13, 14], and shed light into a discrepancy between the two. 
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Figure 1: An illustration of the connections between hydrodynamics, memory matrices, and holography. 

Furthermore, our results can, in principle, be applied directly to microscopic models appropriate for the 
cuprates [6, 7, 15, 16], as we will describe more completely in Section 7. 

1.1. Transport Without Quasiparticles 

Let us briefly review our current understanding of theoretical frameworks for describing transport in 
strongly correlated systems without quasiparticles in two (or more) spatial dimensions. An overview is 
also provided in Fig. 1. 

The simplest way of describing transport without reference to quasiparticles is, in fact, a very old 
framework - hydrodynamics. The modern understanding of hydrodynamics is that it describes the long 
wavelength, long time dynamics of a system close to thermal equilibrium, when there are a small number 
of conserved quantities [17]. We will be liberal, and also allow for some of these conserved quantities to 
decay on long time scales, while still maintaining the name “hydrodynamics”. All we insist upon is that 
this list of conserved quantities is finite; this is in contrast to the typical condensed matter paradigm of 
a Fermi liquid, where occupation numbers at every single wave number are long lived quantities. For 
non-Fermi liquids, it is believed that generic higher dimensional theories do not admit infinite families of 
(nearly) conserved quantities at strong coupling. 

Hydrodynamics proves to be a very powerful framework for describing the dynamics of systems without 
reference to quasiparticles. However, we must emphasize that hydrodynamics is an incomplete description. 
It provides a set of constraints that any (to date) reasonable (2-|-l)-dimensional quantum field theory at 
finite density and temperature, which is approximately translation invariant, must obey, and a universal 
framework within which we can interpret a microscopic calculation. But it does not give us particular 
values or temperature dependence for any microscopic coefficients. 

So now we need a way of obtaining microscopic coefficients. A traditional framework for doing this 
is the memory matrix framework [9, 10]; see [11] for a thorough review. We will review it in this paper, 
but for now we emphasize that it is in principle an exact, microscopic calculation. Its usefulness is that 
it can be efficiently approximated in a hydrodynamic regime where there are only a few quantities which 
do not quickly relax to thermal equilibrium. 

The main purpose of this paper is to clarify and sharpen the connections between the memory matrix 
framework for transport, and hydrodynamic descriptions of transport. In particular, we will not need 
to add a phenomenological momentum relaxation time in hydrodynamics, and compute this coefficient 
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separately using the memory matrix formalism. We will also point out the microscopic computations 
which allow us to compute all phenomenological coefficients within hydrodynamics. 

There has also been, in recent years, a third framework which allows for transport computations: 
gauge-gravity duality. This allows one to study a strongly coupled “large N matrix” quantum field theory, 
directly at finite temperature and density, in real time, by studying a classical gravity theory in one higher 
dimension, in particular black hole backgrounds. Ref. [18] has recently shown the equivalence between 
holographic and memory function calculations in zero magnetic field, in the regime where momentum 
relaxation is slow (this provides the bottom arrow in Fig. 1). One advantage of holographic methods 
is that the calculations are formally valid beyond the regimes of validity of hydrodynamics or memory 
functions. However, we point out in this paper that a wide variety of holographic results for transport, 
computed in regimes where hydrodynamics need not be valid, can nonetheless be understood within the 
framework we derive. 

1.2. Outline of the Paper 

We will begin in Section 2 by recalling the main hydrodynamic results of [5] for electrical transport in a 
simple manner, and using notation suitable for our memory matrix approach. The basic memory matrix 
formalism will be introduced in Section 3, along with its physical interpretation. Section 4 contains our 
main new results on magnetotransport in the memory matrix formalism under the simplifying assumption 
that the diffusion of heat has decoupled from the diffusion of charge and relaxation of momentum: this 
simplification will help us to elucidate many of the subtleties that arise in the memory matrix framework 
when describing magnetotransport. Section 5 relaxes the assumption that heat diffusion decouples, and 
Section 6 describes thermal and thermoelectric transport. 

2. Hydrodynamics 


Let us begin by computing electrical transport in a generic quantum field theory without quasiparticles, 
describable by hydrodynamics, following the approach of [5]. More precisely, we focus on a theory with 
exactly conserved energy and U(l) charge, and approximately conserved momentum, where the only 
dynamics on long time scales are momentum relaxation on the time scale r, and charge/energy diffusion 
on the time scales A^/D, with A the wavelength of the fluctuation and D an eigenvalue of the generalized 
diffusion matrix coupling heat and charge diffusion.^ 

We wish to compute the electrical conductivity matrix, defined as 

Ji — ( 1 ) 

We assume that we have an isotropic, parity-symmetric theory in two spatial dimensions, which constrains 

^xx — 

^xy — ^yx 

Momentum conservation within a hydrodynamic framework - accounting for momentum relaxation - 
gives us the equation 

dtn^ + d^P = -^ + QE, + Beij E (3) 

where Q is the background “charge” density, r is the momentum relaxation time, B is the external 
magnetic field, P is the pressure and H^, is the x-momentum density; a precise definition of Q, applicable 

^Note that with momentum relaxation, there is no longer a propagating sound mode as A: —>■ 0, so indeed the dynamics of 
charge and energy become governed by diffusion. 


(2a) 

(2b) 
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to lattice models, appears later in (27). We have assumed that the velocities relevant for hydrodynamics 
are small and may be treated non-relativistically. As we are interested in perturbing the system with a 
spatially homogeneous electric field, oscillating at frequency w, momentum conservation simplihes to 


iw + 



M5vi 


Q.5Ei + BeijdJj 


(4) 


where A4 is the analogue of “mass density” (namely, the momentum density at small velocities is Ai5vi). 
The memory matrix dehnition of Ai is in (26). The momentum relaxation time, r, is - like all thermody¬ 
namic coefficients - undetermined within hydrodynamics. Unlike thermodynamic coefficients like Ai and 
Q, however, r is extrinsic and sensitive to the precise mechanism by which momentum relaxes, as we will 
see explicitly within the memory matrix formalism. 

(4) gives us an equation relating 6vi, SJi and 6Ei - our perturbatively small quantities within linear 
response. 5Ji may be expressed in terms of 5vi and 6Ei, in an isotropic and parity-symmetric (up to 
magnetic fields) theory, as 

SJi = QSvi + cjq {SEi + BeijSvj). (5) 

Here cjq is an intrinsic “quantum critical” conductivity which measures charge transport independent of 
the momentum mode; a precise dehnition relating it to the charge diffusivity appears later. Combining 
(4) and (5) we may relate SJi to SEi, and we hnd 


_ (r ^ - iuj)AicrQ + .. f ^ \ 

Q^B'^ + ((r“i — iuj)Ai -|- B^aq)'^ \t ) ’ 
_ 2{t-^- iu)AtaQ ++ B^al 
- Q2B2 + ((^-1 _ iuj)At + B^aqf 


(6a) 

(6b) 


The answer simplihes when H = 0 to the form axy = 0 (by parity symmetry) and 


(^xx — T 


Q^r 

Ai{l — iwr) 


(7) 


Though in a system without Galilean invariance one might be concerned that the Qs in (4) and (5) 
need not be the same, we will carefully dehne Q via the memory matrix framework and show that (6) is 
indeed (at leading order) correct. 

These results were hrst derived in [5], for the case of a quantum critical theory, and were argued in 
[8] to be compatible with Hall angle measurements in the strange metal of the cuprates. In fact, the 
structure here is identical, and is valid for any isotropic theory, whether or not it is Lorentz invariant. 
We even need not assume that diffusive charge and heat fluctuations are decoupled (as would be the 
case when the underlying theory has charge conjugation symmetry) - heat transport could only alter the 
equations above if temperature gradients were non-vanishing, which (by assumption) does not happen. 
It is a simple task to compute additional thermal and thermoelectric coefficients, but we will not do so 
using hydrodynamics. 

We remarked previously that the dynamics which governs cjjj is the only dynamics in the system 
occuring on long time scales. However, the expressions for aij are rather complex! This is because there 
are four slow time scales in the problem: the slow driving time scale l/w, and the three time scales 
associated with the specific quantum field theory, and its realization: r, l/wc, and I/wq, where ojc is the 
cyclotron frequency 


Wc 


QB 
M ’ 


( 8 ) 
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and ujq is a “quantum frequency” given by 


Maq' 

Within hydrodynamics, and the memory function formalism, we only work to leading order in a; - but to 
all orders in the dimensionless numbers ujt, ujjLUc, oj/ojq. 

As scales with the strength of the coupling to the operator which breaks translational symmetry 
(as we will explicitly see later), and Wc ~ -B, it is evident in what physical limit these two quantities would 
be perturbatively small. It is less clear how ujq would generically be perturbatively small. One mechanism 
is that the densities of the theory (Q and M) are small (with QjM., the strong-coupling analogoue of 
the fixed charge-to-mass ratio of the electron, fixed). Another is that cjq is anomalously large - though 
this should not be a consequence of weak coupling (i.e., long-lived quasiparticles), for our framework to 
be valid. 

We should also comment that recently, this result has been obtained holographically at w = 0 [8] 
(though it is not expressed in terms of thermo/hydrodynamic quantities) in a set-up which is formally 
valid even when r is comparable to Tm - the time scale associated with microscopic dynamics, and the 
breakdown of hydrodynamics. The limiting case (7), also at cu = 0, is found in various holographic set-ups 
in [19, 20, 21]. It is an interesting question whether or not this is a pathological feature of the holographic 
mechanism of strong momentum relaxation, or a signature that the “hydrodynamic” regime of strongly 
coupled transport persists far beyond our first expectations. 

3. Memory Function Formalism 


Let us now briefly review the memory function framework. As we emphasized in the introduction, the 
method works best when quasiparticles are not long-lived, and the only conserved (or approximately 
conserved) quantities are charge, energy and momentum. In recent years this has become a method of 
choice for studying transport in condensed matter systems without quasiparticles [5, 22, 6, 7], including 
using gauge-gravity duality [23, 24, 25, 26, 27]. We work in units with h= 1. 

Let us consider the set of operators A, B,C,... in a time-translation invariant theory, and correlation 
functions of the form 

IjT 

CAB{t - t') = (A(t)|B(0)) j d\ (A(t)tB(iA)) (10) 

0 

with averages over thermal and quantum fluctuations denoted in (•••). The Laplace transform of this 
expression can be shown to be 

Cab{z) = w [G^q{z) - G5s(i0)] . (11) 

The retarded Green’s function (in real space) is defined as: 

t) = i0(t)([A(x, 0)]) (12) 

where 0 is the Heaviside step function. The momentum dependence of the Green’s functions have been 
suppressed in (11) as we will only be interested in the conductivities, evaluated at zero momentum. As 
standard, we take z to lie in the upper half of the complex plane. We are only computing thermoelectric 
transport coefficients, for which —)• 0) ~ lOjaAB-, with uab strictly finite; thus for us, ^^^(iO) = 0 

and can be neglected. Up to the overall factor of temperature, we recognize Cab{z) in (11) as a generalized 
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conductivity between the operators A and B, ctab- Some formal manipulations on a Hilbert space of 
operators using the inner product (10) give 


o-ab(z) = i^Cab{z) = xac[M{z) + N - izxjcnXDB, 


(13) 


where XAB is the static susceptibility between the operators A and BA 


1 


XAB — Gab{^ — 0) — —{A\B), 


the memory matrix Mab is defined as 


Mab{z) = 


B 


L = [H, o] is the Liouvillian operator with H the Hamiltonian, q is the projection operator 


(14) 


(15) 




AB 


(16) 


and 

^ab = Xab = -XAB- (17) 

Note that Nab vanishes identically in a time-reversal invariant theory. 

We will show that for the magnetotransport problem of interest in this paper, the memory matrix 
effectively truncates to a 6 x 6 matrix, given as m = M -|- N" — iwy in (56). 


3.1. Hydrodynamic Interpretation 

The usefulness of the memory function formalism arises when the inhnite dimensional matrix M has a 
finite number of parametrically small eigenvalues. This allows us to truncate the matrix to a small, finite 
dimensional object. 

To give some physical intuition into (13), let us consider the dynamics of long-lived conserved quantities 
from the viewpoint of hydrodynamics. Let Xa denote quantities which are conserved on very long time 
scales, and let Ua be the conjugate thermodynamic variable. For simplicity, we define Xa by shifting Xa 
by a constant so that in equilibrium, {Xa) = 0. In thermodynamics, we can then relate 

Na = XabUb, (18) 

where y is the same static susceptibility matrix as before, and we drop expectation values within hydro¬ 
dynamics. The transport equations read 

Xa =-MabUb - NabUb + Fa, (19) 

where Nab, as before, is the static susceptibility between A and B, Mab is a (for now arbitrary, but 
symmetric) matrix which denotes the weak relaxation of Xa into an external bath, and Fa is an external 
driving force. Now, we relate Fa to the “thermodynamic field” Ea which drives Xa by 

Fa = XabEb- (20) 

^Differing conventions appear in the literature [28]. This dehnition of the (overall sign of the) retarded Green’s function 
is much more standard [24, 26] and is the common choice in condensed matter physics. 
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Ea should be interpreted as “the net Ua” imposed on the system. Combining these equations together 
at frequency w, we obtain 

Xa = o-abEb ( 21 ) 

where aAB is given by (13). 

One aspect of the memory function framework that is now much clearer is the physical role of the 
memory matrix. Defining the memory matrix as 


Mab = 

in an undriven, time reversal invariant theory (where = 0), we find 

— Tj^j^Xb- 


( 22 ) 


(23) 


Evidently, the memory matrix encodes relaxation times for the average values of the operators A of 
interest. It is now evident why we choose the operators to be almost conserved quantities. In 

this case, the matrix will have parametrically small eigenvalues, and as the memory matrix enters 
the conductivity via an inverse, these quantities will dominate the conductivity. 

This derivation is rather abstract, so let us give a simple example to clarify the thermodynamic 
structure of these equations: the dynamics of a point particle of mass m on a one dimensional line, 
subject to weak friction, but no potential energy. In this case, the momentum p is long lived. The 
thermodynamic conjugate variable is the velocity v, and so the susceptibility is simply the mass. Npp = 0, 
as there is no generic proportionality coefficient between p and v, and so we are left with p = F — MppV. 
Of course if Mpp = 0 we simply have Newton’s Law, F = p = mv\ otherwise, Mpp ^ 0 provides a frictional 
force which dissipates momentum into an external bath. 

It is also important to understand the role played by F in this equation. To compute the “momentum 
conductivity” Opp of this particle, we can interpret the external driving by the force F as actually external 
“shaking” of the line with acceleration Oext = F/m. The conductivity relates the response of p to an 
external source of the time derivative of its conjugate variable. We find at zero frequency: 


p 


mv 


mF 


m 


®ext ®ext Mppttext XIpp 


(24) 


in accordance with (13). 

Of course, we are probably more familiar with electrical conductivity. In this case, the thermodynamic 
conjugate variable to J is a gauge potential —A. And as —dA/dt = E, the electrical conductivity ajj 
relates the response of J to an external electric field E. In this case, as we will see in the next section, 
the analogue of (19) become the hydrodynamic equations which describe the dynamics of charge and 
momentum, which are a bit more cumbersome. 

The usefulness of the memory function formalism beyond the hydrodynamic reasoning described above 
is that we now have explicit expressions for the matrices M, N and y. In particular, this will allow us 
to explicitly compute r, the momentum relaxation time. We will describe more in the following section 
what the components of this matrix actually look like. 

4. Dynamics of Charge and Momentum 


Let us now discuss how the memory function framework produces the magnetohydrodynamic electrical 
transport described earlier. For simplicity, we assume that charge and heat diffusion are decoupled 
processes at zero charge density - we relax this assumption in Section 5. There are two slow time scales 
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in the problem: the time scale over which momentum relaxes tmom ~ t, and the time scale over which 
charge fluctuations diffuse away, tdiff ~ ID, with A the wavelength of the fluctuation; on the longest 
length scales, tdifr oo. Note that at finite charge density, there will be heat diffusion; however, the slow 
diffusion of heat will not be linearly independent from the dynamics of charge or momentum. Our goal 
is to compute the electrical conductivity matrix aij at zero momentum, and small but finite frequency 
w. There are four operators associated with the long lived quantities above which have spin 1 under 
the SO(2) isotropy of the system: the momentum vector P, and charge density fluctuations Vn. This 
suggests that our memory matrix should be a 4 x 4 matrix, as we should only keep track of the (zero 
wave vector components) of these operators. Using the notation of the previous section, our goal is to 
compute the conductivities crj^j^ = axx and aj^jy = axy 

4.1. Diffusive Transport 

Let us begin with the case when the charge density Q = 0, and the magnetic held i? = 0. As i? = 0, 
the time reversal symmetry breaking matrix N will vanish. In this case, hydrodynamics tells us that the 
conductivity is governed solely by diffusive processes. Let us see how this arises in the memory function 
framework. In this case, using the isotropy of the system, in order to compute (Jxx, we need only keep 
track of operators with long time dynamics which are spin 1, and oriented in the x direction: these are 
dxTi and Px - both evaluated at zero wave vector. 

To get started, let us compute the susceptibility matrices. Formally, we introduce a velocity v, as a 
source term conjugate to the conserved momentum, so that Hamiltonian maps as 

Ho^Ho-VxPx, (25) 

where Hq is translationally invariant, and thus commutes with Px, the total momentum in the x-direction 
(we will add terms responsible for momentum decay later). We may compute the susceptibilities XAP,„ by 
studying {A)^^ - the response of the average value of A to a small external source of velocity [17]: 

(26) 

where V is the spatial volume of the sample. This formally defines the value of AI; alternatively, Ai is 
related to Gp p via (14). Note that in the earlier relativistic formalism, Ai was equal to the sum of the 
energy density and pressure, up to coefficients proportional to the (effective) speed of light. Analogously, 

_ {Jx)v,, _ p, 

\PxJx Tr 

Dx ^ 

defines the charge density Q. As the perturbation (25) does not break translation invariance in n(x): 

= 0. (28) 
Identical results hold for the yy susceptibility matrices. By spatial parity symmetry, when H = 0, we find 

XPxPy ~ Xd^ndyTi — XPxdyU — XPydxn — 0- (29) 

In fact, we now see there is a great simplification when Q = 0, because the P indices of the memory 
matrix can play no role in the computation of axx, as XPxJx ~ 9- This is a bit fast; we will see explicitly 
later that the necessary components of the memory matrix Mp^a do not vanish fast enough as Q —)• 0 to 
spoil this argument. 

Next, let us compute XJxdxu- The simplest way to proceed here is to choose a finite momentum regula¬ 
tor, and then show that the transport coefficients are independent of the choice of regulator, kx- We choose 



{Px)vx — \ A 

= Ayr = 
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n(x) = —n(kx) sin(kxx),^ and this implies that the overlap with the current operator Jx(kx) cos(kxx), is 
given by 

(J) = -crqV/r, (30) 

as the conjugate field to n is the chemical potential /u (in the presence of long-range Coulomb interactions, 
should be replaced by the electrochemical potential [29, 30]). We conclude 

XJxdx-n = CTqkx- (31) 


Crudely speaking, we think of fixing x and sending —)• 0 - one then finds that n{x) —)• —n{kx)kxX and 
Jx Jx{kx), so that we have a linear density fluctuation sourcing a constant, diffusive current. But it is 
easiest to compute the memory matrix with the regulator explicitly finite. 

It is a standard result that [17] 


G'nn(k, w) 


Dk"^ — iuj ’ 


(32) 


for any quantum field theory which has a hydrodynamic limit, where n will obey a simple diffusion 
equation at long wavelengths with diffusion constant D: dtn = D'V'^n (recall we have assumed charge 
and heat fluctuations decouple at Q = B = 0)] for the case with long-range Coulomb interactions, (32) 
refers to the irreducible density correlator i.e. the polarizability [29, 30]. We may now use (13) to obtain 
By time reversal symmetry when B = 0, {dxn\dxn) = {dxn\Px) = 0; thus we take the abstract 
projection operator q = 1 in (13). Using that [11] 


Mab{z) 



and regulating this integral, we obtain 


^dxndxni^) 


Dk^ — \uj 


^Q^X + '^^Xnn + • • • ) 


(33) 


(34) 


where we have expanded to lowest non-trivial order in ui, and 


Xnn = ^ = lim lim Re [^^^^(k,^)] 


(35) 


is the static density-density susceptibility. The memory matrix formalism now gives 


^.T.T - 


^Jx dx 77 


^dxTidxn i^Xnn 


Q a; I /n 

= Oci + 0[uj ) . 


(36) 


This is exactly what we should expect ~ the diffusive hydrodynamic coefficient Uq is well-known to be the 
electrical conductivity at Q = B = t). And importantly, the regulator kx has decoupled from the answer. 
We may thus safely take the limit kx —)■ 0. By spatial isotropy, XJxdxn = XJydyn, etc., and XJxdyn = 0, etc. 


4.2. Momentum Relaxation 

Next, let us consider the case where Q ^ 0, and B = 0. We have already computed all relevant static 
susceptibilities, so only the memory matrix M remains. Spatial isotropy implies that Mp^p^ = Mp^g^n = 
0, etc., and so we only need to compute Mg^nPx Mp^p^. Note that the diffusive form of Mg^ndxu and 

®The choice of sine waves instead of plane waves keeps things manifestly real, and so is a bit simpler. 
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XJxdxn from the previous subsection remains valid even with Q 7 ^ 0 , as momentum relaxation kills sound 
propagation at the longest wavelengths, and charge is still conserved, though the value of cjq may change. 

We begin with Mp^p^. Let us consider for simplicity the case where the Hamiltonian of the system is 
given in d spatial dimensions by 

H = Ho- j d^x h{x)0{x), (37) 

where Hq is a translationally invariant Hamiltonian, and 0{x) is an arbitrary operator in the quantum 
field theory. h{x) is an arbitrary function which should depend on space, and be non-trivial in (almost) 
all of the plane. We assume that h{x) is small, so that translational symmetry breaking is weak - we 
will see that this is equivalent to the momentum relaxation time r being large. The operator Pi is easily 
computed to be 

Px = \[H,Px\ = - J d'^x h{x){diO){x) (38) 

The generalization to the case where multiple operators couple to x-dependent fields is straightforward. 
We find 


Mp^Pxiz) 



'■"(gfep.W) 

ui{uj — z) 



Im (Ggg,(ki,k2,a;)) 
OJpJ — z) 


h{]^l)h{^^2)klxk2x■ 


— 00 


(39) 


As before, we have set q = 1 in this equation: we assume the operator O is independent of dxU or Px- As h 
is a small parameter, we may safely take to be the Green’s functions associated with Hq; translation 
invariance then implies that it vanishes unless ki + k 2 = 0. One hnds [5, 24, 27, 18] 


MPxPxid) 


M 

T 


lim 


d^q\Kq)\% 


2..2lm (g0o(q,^)) 


LO 


(40) 


up to a re-scaling of /i(q) to absorb a factor of the spatial volume [18]. We will shortly see that the 
T defined here is equal to the r defined within hydrodynamics; with the memory function framework, 
however, we now have an explicit expression for r in terms of the fields h. A similar formula holds for 
MpyPy. For simplicity, we will assume that the system is isotropic and so Mp^p^ = Mp^p^. 

We have focused on the cj = 0 limit of Mp^p^{lo). This can be justihed on very general grounds. 
Suppose that is the microscopic time scale associated with the quantum dynamics of the system - for 
example, in many quantum critical theories, we have Tm = l/T. The Green’s function Gqq is associated 
with the quantum dynamics of the Hamiltonian Hq, and thus we obtain 

Gooi‘^) = ^00 + (‘^'''m)^oo + (<^T‘m)^^oo -!-••• (41) 

where, for a generic quantum field theory, the coefficients etc., have no anomalously large coefficients, 
relative to any others. This implies that 

M 

^kPxPx {^) = ~ [1 + CujTm H-] , (42) 

with C an 0(1) constant, and thus the hnite frequency corrections to Mp^p^ are higher order in pertur¬ 
bation theory. At cuTm ~ 1, hydrodynamics and the memory function formalism (when truncated to a 
finite set of operators) both cease to be good approximations for generic theories. 
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Next, let us discuss In the case where the operator O is charge conjugation symmetric, this 

matrix element vanishes by charge conjugation symmetry. We discuss the more general case in Appendix 

A. 

Putting everything together, we obtain axy = 0 (by parity symmetry) and, since the memory matrix 
is diagonal up to this point: 


^xx — 


^Jx dx 


+ 




Mg^nd^n - iwXnn Mp^p^ - IxXP^P, 


— (Tq + 


qv 


M .{1 — \ujt) 


(43) 


In a theory without momentum relaxation, we have r —>■ oo; then the real part of the second term in 
(43) is proportional to 5{uj). In such a theory we define cjq to be equal to axx minus the delta function 
contribution, and this definition can be used to compute cjq in a particular model of a strange metal. 


4.3. Magnetic Fields 

Finally, we allow for B ^ We still demand that B is perturbatively small, so that Wc is a perturbative 
parameter. As each component of the memory matrix is first order within perturbation theory (we 
will see this does not change when B 7 ^ 0), it will suffice to consider only the ^-dependent corrections 
to the memory matrix - considering i?-dependent corrections to the static susceptibilities is a higher 
order correction. Furthermore, within the memory matrix, any i?-dependent correction to a parity even 
coefficient such as Mp^p^ must be 0{B^), which is second order in perturbation theory. The only matrix 
which may admit first order corrections within perturbation theory that are linear in B is the time-reversal 
non-invariant matrix N. 

Let us consider the consequences of i? 7 ^ 0 on the matrix N, which relates to static susceptibilities. 
The consequences of an external magnetic field are that 

P, = 5e,,Jf* + --- , (44) 

with the spatially integrated momentum current, where • • • includes effects such as momentum 
relaxation, that we have previously accounted for. The expression (44) represents the Lorentz force law, 
and follows similarly to (38), though with some subtleties. The (zero momentum component of the) 
canonical momentum operator 'Pi(x) which generates translations is no longer equivalent to the physical 
momentum: 

Vi = Pi + j d^x nAp^i, (45) 

where is the gauge potential due to the externally imposed magnetic field. Note that only Pi is gauge 
invariant. In addition, the Hamiltonian is altered to P[ ^ H + Hp, with 


Hp = - d^xJ-Ap. 


(46) 


Now we evaluate Px using the convenient gauge choice Ap = —Byx: 


Px — i[-Po + Ppi Px\ — i 


Hq + Hp,Vx + / d^x nBy 


= \[Hp,Vx]+ / d^xiiBy + 0{B^) 


= - I d^x {-dxJx -h)By = - / d^’x dyJyBy = / d^x BJy = BJy. 


(47) 


As P is a constant, this relates the total current in the y direction to Px, as claimed in (44). A similar 
argument works for Py, with the gauge Ap = Bxy. 
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We can use (44) to find that the B dependent corrections to the matrix N, as xp a ~ ^XJyA for any 
operator A, as an example. Using this logic, we find 


^P:.Py = -^PyP^ = Xp^p^ = -Bxj^P^ = -BQ, (48a) 

^dxTlPy ^dyUPx BXdyUjy BkxCQ^- (48b) 


In writing down these results we have used that XP^d^n = XPydyn = 0, which follows from translation 
invariance of the thermodynamic state. 

To recap: the “full memory matrix” m(a;) at this order in perturbation theory is 


m(a;) = M{u)) + — icox 


/ klaa 

0 

0 

-BkxCTQ 

0 

k^aQ 

BkxCTQ 

A4 ( -iw ) 

0 

0 

-Bkxaq 

-BQ 



A4 ^- iuj 

^ BkxCTQ 

0 

BQ 


\ 


/ 


(49) 


Along with xj^Px = XJyPy = Q, XJxOxn = XJydyu = cfqkx, being the only non-vanishing static susceptibil¬ 
ities, it is straightforward to invert m and obtain the main results displayed in (6), using (13). 


5. Dynamics of Charge, Heat and Momentum 


Let us now relax the assumption that heat transport has decoupled from the problem, and compute the 
electrical conductivity using the memory matrix formalism. In fact, regardless of how many other long 
lived scalar quantities we have - governed by diffusive transport - we will see that (6) continues to hold.^ 
In practice, the only additional conserved quantity is heat. 

Let us suppose that we have currents J/? associated with long lived densities nij(x); the chemical 
potentials associated to these long lived densities are pp] the indices RS ■ ■ ■ will refer to the conserved 
charges, and indices ij ■ ■ ■ will refer to spatial indices as usual. Analogous to (30) we find 

Jr = -'BpsXJpsj (50) 

with S a conductivity matrix. One finds analogously a susceptibility matrix XRS^ and a diffusion matrix 
Drsi with J_R = —Dps^P'S, obeying the Einstein relation 

S = Dx. (51) 

Here we are multiplying together matrices with RS indices. These results are reviewed in [31]. The 
analogous Green’s function to is 

= {klD - iw)^^ klJ^TS, (52) 

and so an analogous computation to before gives us 

^dxnRdxns ~ k^^PS (53) 

Let the index Q denotes the conserved charge; thus 

= Sqq. (54) 

^Indeed, such conserved quantities cannot at zero momentum have any overlap with the electrical current, as the latter 
is a vector and we have assumed rotational invariance in this paper. 
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Denote with P the projection operator onto the Q index: i.e., 

cjqP = PSP. 


( 55 ) 


m generalizes to 




m 


0 

0 




x^QR 


0 

klT.Rs 

-Bkx'BqR 

0 


0 

Bkx^RQ 


-BkxT^RQ \ 
0 


M 


1 


— lo; 


r 

BQ 


M 


-BQ 

1 . 


— lUJ 


(56) 


Our goal is to compute 


-1 


(57) 


and to simply prove that this matrix product is independent of any additional conserved quantities. (It 
is straightforward to carry out the remainder of the algebra and explicitly recover the conductivities.) 

A rather generic and useful result for us will be following formula from linear algebra. Suppose that 
we have a block matrix 

t/=(z w)’ 

with A an m X m matrix, Y an m x n matrix, Z an n x m matrix, and W an n x n matrix. Then 

{X -YW-^Z)-^ -X-^Y{W - ZX-^Y)-^ 


[7-1 = 


-W-^Z{X -YW-^Z)-^ {W-ZX-^Y)-^ 

{X -YW-^Z)-^ -{X -YW-^Z)-^YW-^ 

-{W - ZX-^Y)-^ZX-^ {W-ZX-^Y)-^ 


(59) 


A tedious set of calculations, repeatedly imploying the identities above, leads to the following block 
matrices in m~^: 


mp.p^ = N {B'^ctq + M{t ^ - io;)) 6ij + AfBQcij, (60a) 

(60b) 

^dlnnPj — (60c) 

^ - VB 2 (x(t-i - ii^) + B 2 „o)Phs) Si, - MB^gtRse,,] . ( 60 d) 

where 

N= -^-.. (61) 

Q^B‘^ + {B‘^(j^ + M{T-^-\oj)f 


To obtain (60d), the following identity, along with (55), is helpful (below c is any constant): 

(S + cSPS)“^ = S“^--P. (62) 

1 + CCJq 

Using (50), we find the electric current - conserved density susceptibilities to be 

XJidjTin — ^ij^QR^x- (63) 

It is now a simple matter to see from (57) that the only component of Tirs that enters the final answer 
is Yqq = Uq; thus the magnetohydrodynamic result for the electrical conductivity matrix, given in (6), 
is unaltered by the presence of additional conserved scalar charges. 
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6. Thermal and Thermoelectric Transport 


It is now clear how to extend our work to study thermoelectric transport. When describing thermoelectric 
transport, one must specify a set of three matrices: aij from before, the Seebeck coefficient(s) Oij, and 
Rij, defined by: 


Ji 

Q^ 


^ij 

Ta. 


a. 






-djT 


(64) 


with Qi the heat current. We follow the notation of [5], where R is defined as the linear response coefficient 
between the heat current and the temperature gradient at vanishing electric field; often in experiments 
one measures the thermal conductivity Kij, the coefficient between the heat current and temperature 
gradient at vanishing electric current, J* = 0. The two are related via 


K = K — Taa ^a. (65) 

In the memory matrix framework, we have: 

1 (66a) 

^ijE = XQ.A^ABXBQj • (66b) 

The new susceptibilities we need relate to the heat current. Letting the index H (in the notation of 

Section 5) denote the diffusive scalar quantity heat, we obtain 


XQidjtin — 


(67) 


There is also a new susceptibility: 

XQiPj = 6ijTS. 

This serves as a definition of the new quantity, S, which plays a role of “entropy density”, 
results of (60) thermoelectric transport coefficients may now be computed. Denoting 


( 68 ) 


Using the 


^HQ — chqT, (69a) 

= RqT, (69b) 


we obtain 




^xy — ^ 

f^xx — ^ 


_ {M{t ^-iuj) + B‘^aci)aQ + SQ .. f ^ ■ \ 

_ otQQM{T~^ - iw) + S{Q^ + + aQM(T~^ - ii^))^ 

Q'^B'^ + {B‘^aQ + M{t~^ — iuj)Y 
__ (g^o-Q + A4(r~^ - iu;))(5^ - B^aD - 2SQac:,B‘^ 
Q^B^ + iB‘^a^ + M{T-^-iuj)Y 


QS^ - B^Qal + 2Sa^{B^a^ + M{T-^ - iu;)) ^^ 

Q^B‘^ + {B^aci+MiT-^-ico)Y 


(70a) 

(70b) 

(70c) 

(70d) 


Together with (6), this forms the main result of this paper, and is the most general framework for 
magnetotransport to date, applicable both to relativistic and non-relativistic systems. 

If we specialize to a Lorentz-invariant quantum critical system deformed by a chemical potential (for 
charge) /r, we may compare to the results found in [5] using hydrodynamics. Using (note that we are 
setting the effective speed of light to be 1) 


Ecxq — —^<Tq, 


(71a) 
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TKq= /iVq, 

TS = M- 


(71b) 

(71c) 


we find agreement with the results of [5]. On the other hand, if we set Oq = Kq = 0, then these results 
agree with recent holographic calculations performed at cj = 0 [12, 13, 14] - see also [32, 33]. This sheds 
light into the differences between these papers and hydrodynamics, and suggests how one might go about 
resolving this issue. 

It can be helpful to study Kij instead of Rij ~ in the limit B = 0, Q 7 IO, r—)' 00 , the former is not 
singular while the latter is (see, e.g., [34]). Kij can easily be computed using (65), given our main results 
( 6 ) and (70): 


Kjy 


— ^yy — 4" 


Q S^ (Tq - 25Q'^aq - 
Q-^B^al + (Q2 + M(r-i - iuj)a^y 


rT 


A4(r~^ - icu)(crQ5^ - 2Q5aQO-q - aqQ^) - (TqaqA4^(r~^ - iujf ^ 

Q^B^al + (Q2 + M(r-i - 
(aqQ - crq5)2 


Q^B^al + (Q2 + M(r-i - ia;)aq)' 


^BQT 


(72a) 

(72b) 


Indeed in the r —)• 00 limit, Kxx is finite (only the first line contributes) so long as Q 7 ^ 0. Note that if we 
first set Q = 0 , then Kxx does become singular when r —>■ 00 . 


7. Conclusions 


This paper has developed a model of electric, thermal and thermoelectric transport in strange metals 
which focuses on the influence of a long-lived momentum mode, along with the diffusion of charge and 
heat. Such a long-lived momentum mode is found in essentially all condensed matter models of non-Fermi 
liquids, including those obtained from lattice models appropriate for the cuprates [ 6 , 7]. The influence of an 
external magnetic field on such a mode is universally determined by a few thermodynamic susceptibilities: 
this was established here by the memory matrix formalism, which can be applied to realistic models of 
the cuprates. 

Our results are not valid for systems with spontaneously broken global or gauge symmetries, such as 
superfluids or superconductors. In this case, the Goldstone modes associated with the broken symmetry 
must be consistently included within hydrodynamics, so we expect that they must also be included within 
the memory matrix. Undertaking such task would be an interesting generalization of the present work. 

The resulting B and T dependence of aij, aij and Rij was then reduced to the T dependence of the 
momentum relaxation time r, thermodynamic susceptibilities, and diffusive transport coefficients derived 
via Einstein relations: Uq, aq and Rq. An important feature of this approach is that very different 
physical processes control the values of r and cjq, ccq and Rq . Blake and Donos [ 8 ] argued that reasonable 
assumptions for the T dependence of r and Uq (in particular, r l/r 2 and fjq ~ l/T) lead to an 
appealing explanation of the data on the Hall angle on the cuprates [3]. 

The B dependence of Cij, aij and Rij is explicit, as H is a perturbatively small parameter in this 
framework. Their and our discussions have implicitly ignored the B dependence of Cq, but this could be 
important for understanding experimental thermo-electric data [35, 36]. 

In combination with other recent studies [ 6 , 7, 15, 16], our results now provide a route to the compu¬ 
tation of transport properties of strange metals using microscopically realistic models. 

(i) For the theory of the onset of spin-density wave order in metals, there is a clear separation of the 
degrees freedom responsible for the two terms in (7). The “lukewarm” regions of the Fermi surface far 
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from the “hot spots” contribute to the second term in (7), associated with the slow decay of the momen¬ 
tum mode: a computations of the values of Q, and r in this framework was provided in Ref. [7]. In 
contrast, the “intrinsic” quantum critical conductivity, ctq, in the hrst term of (7) is a property of the 
particle-hole symmetric hot spot theory; the scaling limit of this theory has Q = 0, and so it provides a 
direct computation of cjq [16]. 

{a) The same separation between the lukewarm and hot regions of the Fermi surface applies also to the 
Higgs critical theory of Ref. [15], with their respective contributions leading to the two terms in (7). 

{Hi) For the case of the nematic critical point in two-dimensional metals, the crucial role of the total 
momentum mode was discussed in Ref. [6]. The field theory of this critical point was employed to com¬ 
pute Q, Ai, and r, with weak disorder providing the source for momentum relaxation. There is not yet 
a complete understanding of the value of cjq in such models. 

With these B = 0 computations of Q, At, r, and cjq in hand, then our present analysis shows that the 
extension to weak B ^ 0 in (6) follows immediately, and can be made on quite general grounds. 
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Appendix A. The Memory Matrix Element 


Here we discuss the matrix element when the operator O which couples to the translational 

symmetry breaking held is not charge conjugation symmetric. From the dehnition of the memory matrix. 


M, 


PxdxTL 


= ^ y dw y d'^q h{q) 


Im(-ia;Ggg^„(q, k^,u})) 


UJ[U — z) 


(73) 


where we have related Px to O. Let us discuss the general structure of this Green’s function. In fact, 
noting that z is just above the real axis, we hnd that 


Mp,a, 


n 


d'^'q h{q)Re (Ggg^„(q, kx, 0)) . 


(74) 


We know that Gqq^^ ~ h, as this Green’s function does not obey translation invariance, and thus must 
be proportional to h. Furthermore, so long as charge is an exactly conserved quantity, and this Green’s 
function is analytic in a neighborhood of kx = 0, we conclude that ~ hkx- Thus we conclude that 

^Pxdxn ~ h kx‘ 

Next, let us consider the corrections to the conductivity. For simplicity, we focus on the case where 
B = 0, but similar considerations will hold in the more general case. The full memory matrix m is block 
diagonal, and, considering only xx indices, is 


_ / h^B h‘^kxA \ 

^ K^kxA aqkl ) 

with A and B chosen as functions which are not perturbatively small. We find 

-1 _ 1 f cjq/c^ -h^kxA \ / l/h?B -AjBoc^kx \ 

~ h^klBa^ - klh^A^ V -h^kxA h^B ) \ -AjBa^kx lla^kl J 


(75) 


(76) 
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The correction to the conductivity as given in (7) is 


2AQ 


This is, in our limit, much smaller than either ctq or Ai. Logically, we should think of this as much 
smaller than cjq since we argued previously that either ctq must be anomalously large, or Q must be very 
small, in order for cjq ~ jM. to be possible. (If r is anomalously large and cjq is negligible, then we 
should generically expect boxx to also be negligible). Thus, it is consistent to ignore these corrections to 
the memory matrix at leading order. 
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